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ABSTRACT 

We study the kinematically narrow, low- ionization line emission from a bright, starburst galaxy at 
z = 0.69 using slit spectroscopy obtained with Keck/LRIS. The spectrum reveals strong absorption in 
Mg II and Fe II resonance transitions with Doppler shifts of —200 to —300 km s _1 , indicating a cool 
gas outflow. Emission in Mg II near and redward of systemic velocity, in concert with the observed 
absorption, yields a P Cygni-like line profile similar to those observed in the Lya transition in Lyman 
Break Galaxies. Further, the Mg II emission is spatially resolved, and extends significantly beyond 
the emission from stars and H II regions within the galaxy. Assuming the emission has a simple, 
symmetric surface brightness profile, we find that the gas extends to distances > 7 kpc. We also 
detect several narrow Fe II* fine-structure lines in emission near the systemic velocity, arising from 
energy levels which are radiatively excited directly from the ground state. We suggest that the Mg II 
and Fe II* emission is generated by photon scattering in the observed outflow, and emphasize that 
this emission is a generic prediction of outflows. These observations provide the first direct constraints 
on the minimum spatial extent and morphology of the wind from a distant galaxy. Estimates of these 
parameters are crucial for understanding the impact of outflows in driving galaxy evolution. 

Subject headings: galaxies: ISM — galaxies: halos — galaxies: starburst 



1. INTRODUCTION 

Galactic-scale outflows may be a primary driver of 
galaxy evolution through the removal of cool gas from 
star-forming regions to a galaxy's halo or beyond. Mea- 
surements of outflow properties in distant galaxies, such 
as the mass and energy outflow rates, are central to un- 
derstanding their role in feedback and quenching pro- 
cesses. However, accurate determinations of these out- 
flow rates require knowledge of the spatial extent of the 
wind, a parameter that is difficult to constrain observa- 
tionally. 

While winds are detected in optical line and X-ray 
emission around local sta rburst galaxies to distances of 
sever al kiloparsecs (e.g., iHeckman et all 119901 : iMartinl 
1999), because this emission is faint, outflows in distant 
galaxies have traditionally been detected only in absorp- 
tion against the stellar continuum. This technique probes 
cool ga s (T < 10 4 K) trac ed by transitions such as Na I 
D (e.g.JRupke et al.ll2005l) and Mg II (e.g. JWeiner et all 
l2009tlRubin et aL 201 0q), and constrains the outflow ve- 
locity, covering fract ion, and column density. Recently, 
iMenard et alJ ( 20091 ) also showed that Mg II absorption 
tracks the overall star formation history of the Universe. 
However, these absorption-line studies provide only weak 
constraints on the radial extent, spatially-resolved geom- 
etry, and volume density of the outflow. In turn, esti- 
mates of the rates of mass and energy loss are uncer- 
tain by at least two orders of magnitude. Improved con- 
straints are critical for understanding the role played by 
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outflows in polluting a galaxy's halo and the surrounding 
intergalactic medium. 

A novel technique for studying outflows in the distant 
universe is analysis of emission from the outflowing gas. 
This emission may arise from resonance-line scattering 
off of the flow or fluorescent radiation powered by lumi- 
nous star clusters. In the former case, the emission, with 
the corresponding absorption, may exhibit a P Cygni- 
like line profile. The high optical depth of the Lya tran- 
sition results in P Cygni profiles observed in the spec- 
tra o f Lym an Break Galaxies (LBGs; e.g.. iPettini et all 
I2001L 120021 ). Although it has much lower optical depth 
than Lya, Mg II AA2796, 2803 photons can also be reso- 
nantly trapped and will produce strong emission where 
the gas is optically thick (i.e., at hydrogen column den- 
sities of Nh ~ 10 19 cm" 2 ). In contrast, optically thick 
Fe II resonance absorption is not trapped, but pumps 
UV Fe II* emission lines at ~ 2000 - 3000 A. In princi- 
ple, Mg II and Fe II* emission can be particularly useful 
for tracing winds at z > 0.3, where it is detected into 
the optical; moreover, Mg II emission has already been 
detected in z ~ 1 star-f orming galaxies (|Weiner et al.l 
120091 : iRubin et al]l2010bf ). Measurement of the spatial 
extent of this line emission provides a firm lower limit on 
the radial extent of the outflow and constrains the wind 
morphology. 

In this paper, we examine the line emission from a 
bright starburst at z = 0.694 known to exhibit a P 
Cygni profile in Mg II and Fe II* emission (jRubin et al.l 
l2010a| ). We report in detail on the observed charac- 
teristics of the emission and suggest possible produc- 
tion mechanisms. We adopt a ACDM cosmology with 
H = 70 km s- 1 Mpc" 1 , fi M = 0.3, and Q A = 0.7. 

2. OBSERVATIONS 
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Fig. 1. — The Mg II line profile in the galaxy spectrum. Dot- 
ted portions of the spectrum are associated with the transition in 
the opposite panel. The systemic velocity is marked with verti- 
cal dotted lines, and the horizontal dashed line marks the contin- 
uum level. The profile exhibits blueshiftcd absorption extending 
to —800 km s — 1 relative to the systemic velocity. Emission at and 
redward of systemic velocity is also evident. Together, these fea- 
tures exhibit a P Cygni profile, suggestive of a galactic outflow. We 
propose that the red and blue sections of the spectrum in the left- 
hand panel arise from different areas of the outflow, as indicated 
in Figure [4] 



Details of our observations of the t arget galaxy 
(TKRS4389) and data reduction are given in lRubin et al.l 
(|2010al) . We obtained spectroscopy of this galaxy us- 
ing the Lo w Resolution Imag ing Spectrometer (LRIS) 
on Keck 1 (jCohen et al.lll994j ?. Our instrumental setup 
afforded a FWHM resolution ranging between 190 — 
300 km s^ 1 and wavelength coverage of - 3200-7600 A. 
We used a 0.9" s litlet oriented NE (see Figure 1 of 
I Rubin et al.l [20 10 a) and collected six — 1800 sec expo- 
sures with FWHM — 0.6" seeing. The data were reduced 
using the XIDL LowRedus0 data reduction pipeline. 

3. ANALYSIS 



As noted in iRubin et al.l (|2010aD , this galaxy is excep- 
tionally bright for its redshift, with a star formation rate 
of — 80 Mq yr" 1 . Stellar population modeling indicates 
that the spectrum is dominated by light from the intense 
star formation activity. The deepest parts of the Mg II 
and Fe II resonance absorption lines are blueshiftcd by 
— 200 — 300 km s -1 , with high velocity tails extending to 
—800 km s -1 , indicating that these ions trace an outflow 
(see Figures [T]and [2]). Here we analyze the characteristics 
of the observed emission from these ions. 

3.1. Characteristics of Mg II Emission 

The one-dimensional spectrum of this galaxy reveals 
strong Mg II emission near systemic velocity and extend- 
ing to the red (Figure [J). We measure the kinematic 
extent v exten t of each emission line (Ta ble [1]) using stan- 
dard techniques (jCooksey et al . 1 120081 ) . Line fluxes are 
derived from integrating the continuum-subtracted spec- 
trum over this velocity range. The continuum is mea- 
sured from clean regions near each feature; its error is 
determined through Monte Carlo realizations. The flux- 
weighted velocity centroid of the lines is given as v cen , 
and Vp ea k is the brightest pixel in each line (smoothed 
by 3 pixels; Table 1). 

4 http:/ /www. ucolick.org/~xavier/LowRedux/ 
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Fig. 2. — Fe II transitions in the galaxy spectrum. The sys- 
temic velocity is marked with vertical dotted lines, and the hori- 
zontal dashed line marks the continuum level. The left-hand col- 
umn shows resonance absorption lines, while the right-hand column 
shows Fe II* emission profiles (top four panels) and the [Ne III] 
A3869 line (bottom panel; black). The green line in the bottom 
panel shows the coadd of the detected (and unblended) Fe II* emis- 
sion lines. Gray arrows mark Fe II* transitions arising from states 
that cannot be radiatively excited directly from the ground state 
(see J|4]and Table 1). We propose that the red and blue sections 
of the spectrum in the left-hand panel and the magenta and cyan 
sections in the right-hand panel arise from different areas of the 
outflow, as indicated in Figure [4] 



Mg II absorption and emission are also evident in the 
two-dimensional spectrum of the target (Figure [3]). To 
explore the spatial extent of the Mg II emission, we sub- 
tract a model of the two-dimensional continuum pro- 
file from the original two-dimensional spectrum. The 
residual emission extends up to ~ 2" from the spatial 
center of the galaxy continuum, and extends well be- 
yond the stellar structure (i.e., the components of the 
galaxy which generate the continuum flux near Mg II). 
The flux of the two emission features at 2796 A below 
and above the galaxy continuum are (8.0 ± 0.4) and 
(4.4 ± 0.4) x 10~ 18 ergs cm~ 2 s _1 , respectively, and 
have flux-weighted velocity centroids of — 26 km s" 1 
and 74 km s _1 . The corresponding fluxes at 2803 A are 
(4.0 ±0.3) and (2.5 ±0.4) x 10" 18 ergs cm" 2 s" 1 . There 
is a significant difference in the flux of the 2796 A line in 
these two locations; this may reflect the geometrical dis- 
tribution of the emitting gas, suggesting a lack of spheri- 
cal symmetry. Combined, these features contribute most 
(~ 80%) of the total 2796 A flux. In Figure [3k we also 
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Fig. 3. — (a): Two-dimensional Keck/LRIS spectrum of the Mg II doublet, evident in emission and absorption. The vertical blue and 
red lines mark the systemic velocities of the Mg II doublet transitions, (b): A model of the galaxy continuum, scaled by the flux in the 
spatial center of the object, has been subtracted from the spectrum shown in panel (a). The residual flux at A = 2796 A and A = 2803 
A indicates that the Mg II emission is broader than the continuum, and extends on both sides to ~ 1" (> 7 kpc). (c): The spatial 
profile of the spectrum shown in panel (b ), averaged over the full velocity extent of both emission lines (black circles ). T he lines show the 
continuum-subtracted emission profile from models with different distributions of line-emitting material. See Section l3.1l for details. 



present the continuum-subtracted spatial profile aver- 
aged over the velocity range covered by both lines. Errors 
calculated from the extraction of continuum-subtracted 
spatial profiles at other locations along the slit are con- 
sistent with the error bars in the Figure in the wings of 
the profile and are larger by up to a factor of ~ 2 near the 
profile center. Spatially extended emission from strong 
lines arising in H II regions, e.g., H7 and [O II], is not 
detected. 

Finally, we develop a rudimentary model for the spatial 
distribution of the line-emitting material to constrain its 
distance from the galaxy's center. Because Mg II is emit- 
ted in star-forming recombination regions ()Kinnev et al.l 
1993), the model assumes that the emission is distributed 
over both the galactic H II regions and a more extended 
area with variable geometry and size. The nebular emis- 
sion component is distributed in roughly the same way as 
the young stars in the galaxy; i.e., as in the HST/ACS 
6435-band image. We create a "continuum" model by 
scaling this image such that when it is convolved with a 
two-dimensional Gaussian to simulate the effects of see- 
ing and is "observed" with the correct slit orientation, the 
spatial distribution of the flux passing through the slit 
matches the galaxy continuum profile. We then assume 
that the extended emission is distributed either (1) with 
uniform surface brightness and an outer radius i? ut, or 
(2) in a ring with radius -R r ; ng and thickness ~ 1 kpc. 

We generate several realizations of these models by 
allowing the ratio of the total luminosity from the ex- 
tended emission to the luminosity in the section of the 
galaxy centered on the slit (/ ga i) to vary between 0.5 
and 3 in increments of 0.25, and by varying either i? ut 
or i? r ; ng between 5 and 13 kpc in increments of 0.25 kpc. 
These realizations are added to the 6435 image, and the 
resulting "extended emission image" is convolved with 



a Gaussian and "observed" through the slit. We then 
subtract the profile of the continuum model from the 
extended emission model profile, first scaling the contin- 
uum profile to the peak flux value in the extended emis- 
sion profile. In general, we find that the uniform surface 
brightness model provides acceptable fits to the data if 
1.25 < /gai < 2.5, although adjustments in the value of 
i?out are required. With / ga j = 1.5, i? ut must exceed 9.5 
kpc to match the observations. As / ga i increases, values 
in the range 8.25 kpc < i? ut < 13 kpc provide an excel- 
lent match to the observed emission (Figure |3J;; red solid 
line). Additionally, though the best- fit uniform surface 
brightness model yields a lower x 2 value than any of the 
ring-like emission models, realizations of the latter with 
6.5 kpc < i?ri ng < 8.25 kpc and / ga i ~ 1.25 are allowed 
by the data (gray dashed line). While this framework 
is quite simplistic, the models suggest that the material 
giving rise to the Mg II emission is distributed over a 
large area, extending to distances of at least 6.5 kpc. 

3.2. Characteristics of Fe II* Emission 

The Fe II ion is notable for its multitude of permitted 
transitions in the wavelength range 2000 - 3000 A (Ta- 
bled]). We observe several Fe II resonance lines in absorp- 
tion (Figure [2]) which are blueshifted by ~ —200 km s _1 , 
similar to Mg II. We also identify several transitions in 
emission (Figure[5]). All of these are consistent with tran- 
sitions arising exclusively from the J = 9/2 or 7/2 upper 
levels (see Table Q] for details). 

The velocity limits of each emission line, v cen , and 
Vpeak are determined as for Mg II. Both v pea k and v cen 
are redward or within 30 km s _1 of systemic in 4 of 5 de- 
tected lines. In an inverse variance- weighted stack of the 
emission lines at 2365 A, 2396.35 A, 2612 A, and 2626 A 
(overplotted in green in Figure [2]), v pea k = +50 km s _1 
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Fig. 4. — Cartoon illustrating resonant scattering in Mg II and 
fluorescence in Fe II. The galaxy is surrounded by an expanding, 
spherically symmetric wind. Mg II and Fe II ions in the section 
of the wind closest to Earth absorb continuum photons at the rest 
frequency of the ions; this produces a blueshifted absorption line 
profile (shown in blue in Figures [T] and |2j . Ions in the back side 
of the wind scatter photons with wavelengths redshifted relative 
to systemic velocity; these photons travel through the wind and 
into our line of sight (shown in red in Figures [T] and [SJ . Because 
there are several energy levels accessible to Fe II, the photons are 
not resonantly trapped, and escape from the wind via emission to 
excited states (indicated in cyan and magenta in Figure [2}. 



and v cen — +8 km s -1 . We note that this differs from the 
veiocity profiles of nebular emission lines such as [Ne III] , 
H5, and H7, which have peak velocities near or blueward 
of systemic and blueward of v cen . The [Ne III] velocity 
profile is shown in Figure [2] as a reference; this line is 
chosen because the profile is not contaminated by broad 
stellar absorption. The velocity profiles of the Fe II* 
emission lines significantly differ both from those of the 
Fe II resonance absorption lines and from the profiles of 
lines tracing the kinematics of the galactic H II regions. 

4. DISCUSSION 

We now suggest a possible production mechanism for 
the observed Mg II and Fe II* emission: photon scat- 
tering in a large-scale galactic outflow. As discussed in 
Section [1] this has been invoked in previous work to ex- 
plain the P Cygni profile of Lya emission in LBGs. Sim- 
ilarly to Lya, Mg II is likely to have high optical depths 
in the ISM of star-forming galaxies. Additionally, the 
Mg II transitions of interest arise from a set of two close 
upper energy levels that can only decay to the ground 
level. Mg II photons are therefore resonantly trapped in 
high optical depth conditions (see also Prochaska et al., 
in prep). 

Figure 0] illustrates the effects of a transition reso- 
nantly trapped in an outflow. Mg II ions in the sec- 
tion of the wind closest to Earth absorb continuum pho- 
tons at the Mg II transition in the rest frame of the gas; 
this produces a blueshifted absorption line profile. Mg II 
ions in the far section of the wind scatter photons with 
wavelengths redshifted relative to the rest frame of the 
front part of the wind; these photons can therefore travel 
through the wind, generating emission at and redward of 
the systemic velocity. Such an outflow may also produce 



the spatially extended emission evident in Figure [3] This 
emission arises from gas above and below the galaxy in 
Figure |4] moving tangentially to our line of sight, and so 
is expected to occur near the systemic velocity. 

A similar mechanism can produce the observed Fe II* 
emission. Ions in the front side of the wind will again pro- 
duce blueshifted absorption in resonance transitions, as 
in the case of Mg II. The excited ions in all parts of the 
wind will then emit photons and decay to states with 
I A J] < 1, including Fe II* transitions to fine-structure 
levels with J > 5/2. We note that these are the only 
transitions observed; i.e., we do not detect emission lines 
with lower energy levels having J = 3/2 or 1/2, and the 
detection of these lines in higher-S /N data would rule out 
this picture. Further, the observed emission arises from 
the same upper levels excited in the resonance transi- 
tions. The lack of observed absorption from these ex- 
cited states indicates that they decay before the ion is 
re-excited. The wind is therefore completely optically 
thin to the observed Fe II* photons. 

5. IMPLICATIONS 

The mechanism described above is only one of sev- 
eral which may produce Mg II and Fe II* emission. 
These lines may be produced in recombination regions 
(IKinnev et all 1 19931) or in AGN (|Vestergaard fe Wilkesl 
120011) . However, because the kinematics of all the lines 
discussed are qualitatively different from the kinematics 
of the nebular emission in the spectrum, and because 
our analysis of the spatial distribution of the Mg II emis- 
sion suggests that it arises from distances out to at least 
~ 7 kpc from the galaxy center, we consider photon scat- 
tering the most probable production mechanism for the 
observed emission. As such, our observation represents 
the first detection of a galactic outflow in emission in 
a distant galaxy, and provides the first constraints on 
the minimum spatial extent f> 7 kpc) and morphology 
of this outflow. Stringent constraints on the mass and 
energetics of the gas require a solution to the radiative 
transfer equation in the wind. However, here we use our 
EW measurements (Table [l) to make a crude estimate of 
the physical properties of the outflow. The ratio of the 
Mg II absorption line EWs (i.e., the "doublet ratio") is 
~ 1.64; this indicates that the optical depth is T273fj > 1.4 
along the line of sight to the galax y itself (|Spitzerlll968l : 
I Jenkins! [l98l iRubin et al] l2010bl ). This analysis does 
not account for the effects of the emission on the absorp- 
tion EWs; however, because the 2803 A absorption line 
is filled in by emission from both transitions, the EW 
is likely reduced more than the 2796 A absorption EW. 
We therefore consider this estimate a conservative lower 
limit. Assuming a coverin g fract ion equal to 1 and in- 
voking Eq. 2-41 of Spitzcr (1968J), T2796 > 1-4 implies a 
column density of A^Mgii > 10 14 cm~ 2 . 

In the case of an isotropic wind which is fully cov- 
ered by the slit, the sum of the EWs of absorption and 
emission in a given transition equals zero (Prochaska et 
al., in prep.). The sum of the Mg II EWs measured 
in the galaxy spectrum is ~ 0.7 A, indicating that the 
outflow is nearly isotropic, and hence that the sections 
of the wind observed in emission have T2796 > 1-4 (i.e., 
similar to the portions of the wind observed along the 
line of sight). With w w i n d = 300 kms _1 , Au w ind = 
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500 km s , and AR wln d = 7 kpc, and using the equation 

72796 = (ne 2 I m e c) ,f2 796A2796nMgIl(Ai? win d/ AUw ind), 

modified from Eq. 8.45 of lLamers &: C assinclli ( f 999|), we 
find a number density riMgii ~ 7x 10~ 9 cm -3 . Neglecting 
any ionization correction, and assuming solar abundance 
(log Mg/H = -4.42) and a fact or of -1.2 dex dust deple- 
tion (jSavage h Sembachl 119961 ). the corresponding limit 
on riH is 0.003 cm -3 . While this estimate suffers from 
large uncertainties, the limits on both the wind radius 
and riH determined from this analysis provide novel in- 
formation which is crucial to constraining both the mass 
and energy carried by the outflow. 

Beyond the analysis presented here, future studies of 
these emission features in other galaxies will provide crit- 
ical insight into the impact of galactic winds on the redis- 
tribution metals and dust to the o uter halos of galaxie s 
(jMenard et al.ll2010h and the IGM (|Simcoe et all [2002). 
A rest-frame UV spectroscopic survey of 0.3 < z < 1.4 
galaxies (Rubin et al. 2010, in prep.) reveals that these 
emission features are quite common at z ~ 1, and thus 



may be targeted for follow-up observations with IFUs to 
create maps of the outflow morphology. High spectral 
resolution observations of Fe II* emission will probe the 
gas kinematics in detail from the outskirts of the wind to 
deep within the galaxy, both on the far and near side of 
the stellar disk; in addition, the observed flux in the fine 
structure lines can provide an in dependent constraint o n 
the radial extent of the gas (e.g.. iProchaska et al.ll200d) . 
These studies, when performed in concert with radiative 
transfer modeling of the emission features, will probe the 
morphology of outflows at an unprecedented level of de- 
tail in numerous distant galaxies. 
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9/2<-9/2 


2374.46 
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Fe II UV3 


42658.22 


0.00 


7/2^9/2 


2344.21 


Absorption 


2.08 ±0.08 












42658.22 


384.79 


7/2^7/2 


2365.55 


5.90E+07 


-0.17 ±0.06 


0.252 ±0.091 
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43238.59 


384.79 


5/2^7/2 


2333.52 


1.27E+08 




< 0.301 










43238.59 
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5/2-5-5/2 


2349.02 


1.09E+08 




< 0.292 










43238.59 


862.62 


5/2-5-3/2 


2359.83 


5.42E+07 




< 0.284 










43620.96 


862.61 


3/2^3/2 


2338.73 
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Mg II 








2796.35 


Absorption 


2.65 ±0.07 






















-1.69 ±0.05 
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1.62 ± 0.05 






















-1.84 ±0.06 
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Ne V] 

O II] 

Nc III] 
H5 
Ho 



3426.98 
3727.10 
3869.84 
4102.90 
4341.69 



-0.61 ± 0.13 
-66.74 ± 0.25 
-5.74 ±0.13 
-4.84 ±0.15 
-14.42 ± 0.30 



0.466 ±0.103 
50.263 ±0.190 
4.673 ±0.111 
3.591 ±0.117 
7.912 ±0.207 



51 
194 
-26 
-39 
1 



212 
132 
25 
38 
81 



-271—5- 695 
-464-5- 732 
-373-5- 436 
-312-5- 397 
-407-5- 621 



Note. — Rest-frame EW, flux and kinematic measurements from LRIS spectro scopy. Upper and lower energy levels, orbital angular 
momenta (J) and Einstein coefficients are given for Fe II transitions {Morton 2003). Errors are lcr uncertainties with 3<j limits. Blended 
transitions are marked with "b". We check Fe II identifications by comparing A coefficients among different transitions. For optically thin gas, 
the flux ratio of two lines that originate from the same excited state is approximately the ratio of the corresponding A coefficients. E.g., the 
flux ratio of the 2612 A and 2632.11 A transitions (Ehigh = 38660.04 cm -1 ) is ~ 2.8 ± 0.7, i.e. ~ 1.2<r higher than the ratio of A coefficients 
(1.98). 

a Line fluxes are determined by summing the continuum-subtracted spectrum over the velocity range of the line, as described in £13. II 



